Two dispersion curves for an interacting Bose gas in a vessel 
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In the framework of the Gross-Pitaevskii approach, we have considered the interacting Bose gas 
in a one-dimensional bounded domain and have found two dispersion laws for phonons. One law 



coincides with the well-known Bogolyubov one hoj(k) 



2m 



+ qnv{k)^- with q = 1. 



The 



second law is new and is described by the same formula with q = 1/2. The first solution corresponds 
to the single harmonic (as for the cyclic boundaries) , and the second solution is represented by the 
set of harmonics forming a wave packet. Two solutions appear due to two different combinations of 
harmonics in the wave packet. 

PACS numbers: 67.10.Ba, 67.85. De 
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I. INTRODUCTION 

The first models of weakly interacting Bose gas 
(WIBG) were constructed by Bogolyubov [l| and by 
Bogolyubov and Zubarev Q. The dispersion law for 
phonons in WIBG was also reproduced in the Gross- 
Pitaevskii (GP) approach [3|— l6|j . After the obtaining of 
a condensate of dilute gases in magnetic traps, WIBG 
was intensively studied for two decades (see, e.g., reviews 
However, one unusual important property was not 
noticed, to our knowledge. It is related to the influence 
of boundaries on the dispersion law. 

Bogolyubov dispersion law was obtained for cyclic 
boundary conditions (BCs). With regard for the bound- 
aries and the inhomogeneous field of a trap, the disper- 
sion law for small k was deduced in works | 9|4llll . and 
it agrees approximately with experiments |12l Il3| . Be- 
low, we will determine the dispersion law for a homoge- 
neous medium and for all k taking the vessel walls into 
account, but without the field of a trap. Such a task is al- 
most equivalent to Bogolyubov one [l|, [2| . The difference 
consists in the consideration of the boundaries, which re- 
sults in the appearance of one more dispersion branch. 
Below, we will deduce two dispersion laws and will study 
their properties. The difference in the solutions for cyclic 
and noncyclic systems is related to the difference in the 
topologies of these systems. 



II. SOLUTIONS OF THE GROSS-PITAEVSKII 
EQUATION 

For simplicity, we restrict ourselves by the one- 
dimensional case (ID). The well-known theorems [l4]-[lj| 
forbid the existence of a condensate in the ID and 2D 
cases at T > 0. This prohibition is valid for infinite 
systems. However, all systems in the Nature are finite. 



It is sometimes neglected this circumstance and consid- 
ered that the prohibition is true for the finite systems as 
well. But this is not quite right: for the finite systems 
with nonpoint interaction, the condensate is permitted 
in ID at very small but finite T (see Appendix). Even if 
the condensate would not be present, we would consider 
the GP equation in the ID case formally, not identifying 
^(x, t) with the condensate; the effect obtained below 
is related to boundaries, rather than to the dimension. 
Therefore, this effect must be present in the 3D case, 
where the condensate is surely present. 

We now pass to the solution of the problem. Consider 
WIBG in a vessel without external fields. The general 
form of the GP equation is 



ih 



h 2 d 2 ^{x,t) 



dt 2m dx 2 

+ •M)/**fl*->W.'>f' 

with the normalization 

dx'\^(x',t)\ 2 = N, 



(1) 



(2) 



where N is the total number of particles. The GP equa- 
tion can be deduced from the Schrodinger iV-particle 
equation in the mean-field approximation [6[ or from the 
Heisenberg operator equation [!, In both cases, one 
can understand ^(x,t) as the wave function (WF) of an 
inhomogeneous condensate. 

Let the field created by the vessel be as follows: 



U w (x) 



U s >0 x< -L/2, 
-L/2 < x < L/2, 
U s x> L/2. 



(3) 
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For such a potential of walls, the WF of the system 
changes sharply its behavior at x = ±L/2. Therefore, 
the solution for ^(x,t) should be sought separately in- 
side and outside the vessel. Then these solutions must be 
sewed together. Below, we will find the solution inside 
the vessel and consider the limitations imposed by the 
sewing conditions. 



2 



In the ground state, the condensate is described by the 
WF 

V (x, t) = R (x)e lS / h , S = -Et, (4) 

which satisfies the equation 

dW _ h 2 d 2 t> 
dt 2m dx 2 

+ 9 (x,t) [ dx'U(\x-x'\)\^o(x',t)\ 2 . (5) 



The condensate is homogeneous in the whole volume, ex- 
cept for the narrow domains near the walls. Its concen- 
tration n(x) = R 2 {x) = no(x). In the presence of small 
oscillations, we have R{x) = Ro(x) + 5R(x) and 

n = no(x) + ho(x, t), S = — Et + sq(x, t). (6) 

Let R, n and S in (g]), © be real. Then relations ©-© 
yield the equations for small ho(x,t) and so(x,t) (which 
should be also real): 



dn 
dt 



-V [Vs (n + n )] , 



(7) 



dt 



h 2 ^[(no+no) 1 / 2 



1/2 



(Vs ) 2 



2m (no 4 
h 2 6R V 2 R 
2m Rq R + SR 



no) 1 / 2 2m 

dx'fio(x' , t)U(\x — x'(\§) 



To simplify the equations, we neglect the inhomogeneity 
of Rq(x) near the walls, by setting Vno = and Vi?o = 0. 
This is admissible, since Vno and Vi?o ar e considerably 
different from zero only at a distance of the order of the 
mean interatomic one from the wall. Therefore, the in- 
fluence of such an approximation on the dispersion law 
must be small. Moreover, we consider oscillations of the 
density and the phase so weak that their smallness ex- 
ceeds the smallness of the potential. Then we can restrict 
ourselves by the linear approximation. Thus, we obtain 



1 dT s (t) 
T n (t) dt 



= C 2 , 



dn(x) = -^\7 2 s(x), 
m 



(13) 
(14) 



C 2 s(x) = -- V 2 n(.x)+ I dx'n(x')U(\x-x'\). (15) 

4mno 



The solution of Eqs. (1121) and (|T3|) can be written as 

T„(i) = coswi, T s (t) = smut, (16) 

d = C 2 = -u, (17) 

in the real form or 

T n (t) - T s (t) = Ci - -C 2 - iuj (18) 

in the complex form. The real values of ho(x,t) and 
So(x,t) are obtained only for solution p^|) . (fT7|) . which 
will be used in what follows. Then the equations for n{x) 
and s{x) take the form 



ujfi{x) — — V 2 s(x), 



(19) 



-ljs(x) = -— — V 2 n(x)+ / dx'n(x')U(\x-x'\). (20) 

These are two basic equations under study. They are 
rather simple, but it turns out not very easy to solve 
them. 

Let try to seek a solution as a single harmonic: 

h{x) = a 2l e^' x , s(x) = b 2l e lk ^ x , k 2l = 2nl/L. (21) 

Under periodic BCs, the potential can be expanded in a 
Fourier series 

U(\ Xl -x 2 \) = j J2 v(k)e lk{xi - X2) , (22) 

k=2-Kj/L 



dh n 2 

-7T- = V s , 

at m 



(9) 



V 2 n + / dx'n (x',t)U(\x-x'\). (10) 



ds _ 
(If AmriQ 



Since only the standing waves can be stationary states 
in the presence of boundaries, we seek the solutions in 
the form 

h (x,t) = h(x)T n (t), s (x,t) = s(x)T s (t). (11) 

Substituting these functions in Eqs. ([§]) and (fTU)) and 
separating the variables, we obtain 



1 dT n (t) 
T s {t) dt 



Ci, 



(12) 



v{k) = I U(x)e~ lkx dx 



(23) 



where x = x\ — x 2l and j runs all integers. Under periodic 
BCs, if xi,x 2 £ [-L/2, L/2], then x\ — x 2 belongs also 
to [-L/2,L/2], Substituting ([2TJ) and (JUJ) in Eqs. ^ 
and l|20p. we obtain the Bogolyubov dispersion law 



E b (k) = 



1 (h 2 k< 



V 2m 



n u{k) 



h 2 k 2 



in 



(24) 



If we will try to solve similarly the problem with regard 
for boundaries, we will fail. In this case, we have £1—0:2 £ 
[— L,L], and the potential is given by the series 



1 



U(\x 1 -x 2 \) = — ^ u(k) 



(25) 
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(j are integers), which follows from (|22|) at L — >• 2L. 
Integral (|23p is now taken over the interval x € [— £, L]. 
For the short-range potential, the result for v(k) is the 
same as in the case of periodic BCs. Under periodic BCs, 
the integral / dx'h(x')U(\x—x'\) in (|20|) is reduced to one 
harmonic. But, in the presence of boundaries, we have, 
with regard for (l2"5")) . 



v Q t j) inj(x-x')/L 

2L 



dx'n(x')U(\x-x'\) 

L/2 

J dx'a 2l e a ^'/ L J2 

-L/2 j 
_ OvMtet) i2nlx/L , a '2l ST^ ,.( h \ 



farfafa+lwz sinK2?-2j a -l)/2] 

2/-2j 2 -l ' V 7 

where j,j2 run all integers. The integral becomes sepa- 
rated into the sum of the initial 2Z-mode and many odd 
modes. Therefore, the solution should be sought as a 
sum of many harmonics. As it will be seen, there exist 
two solutions. 



where j and I are integers. Using the expansion 

e iq 2 j + ix _ e i7r(2j+l)x/£ 

= 2^ sin fr(2p-2j-l)/2] i7r2pg/L 
~ ir ^ 2p - 2j - 1 



(31) 



(p are integers) , we reduce the last sum in (|2"9"j) to a sum 
of the form ^e jfc2!a: j3 ; . Then we gather and equal the 
i 

coefficients of the functions e lk2lX in (|2"9")l to zero. In this 
way, we obtain the equations 

2 ^ , ,sin[7r(2j + l)/2] 

7T » — • 



E 



2j + l 



sin [ tt(2p - 2j - l)/2] 
2p-2j-l 



p#0 



2p-2j-l 



(32) 



a 2 *(^(fc 2 /)-^ 2 ) 
4^ 2 ^ siii[7r(2Z-2j-l)/2] 

^ 4- ^ (<Z2J+1 } 2Z-2j-l 

j 

ST sin[^(2p-2j-l)/2] 



ft 2 A: 2 

4fe J+ i) = n v(q 2:j+1 )^-^. 



A. First solution. 

To obtain it, we represent the functions h(x) and s(x) 
in the form of expansions in the full collection of "even" 
exponents 

n(x) = a 2 ie tk2 ' x , s(x) = J2 b2ie lk2,x , (27) 



where k 2 i = 2nl/L. Substitute these series in Eqs. (IT9"1) . 
(|20f with regard for (p5|) . We note that the exponents 
e zk 2 ix are i n d e p en d en t basis functions. Therefore, the 
equations for the coefficients a 2 i and b 2 i and the fre- 
quency u> can be obtained, if the sums of coefficients of 
each of the exponent e 2lX are set to zero. Thus, relation 
(H]) yields 



ujm 

o-o = 0, b 2 i^ = —a 2 i. 



nok 2 



(28) 



2/ 



In view of (I2U1) and (|2"5)) , we obtain 



-w&o =^a 2 /e lfe ^— ^(i^(fc 2; )-^ 2 ) (29) 
; n n. /c 2; 



a 2 ;- 



sin [tt(2Z — 2j — l)/2] 
tt(2Z - 2j - 1) : 



E}(k) = 



h 2 k 2 

2m 



n Q v(k) 



h 2 k 2 

2m ' 



(30) 



(34) 



Since 



sin [tt(2Z - 2j - 1) /2] sin [tt(2j3 - 2j - l)/2] 
= cos (77/) cos (7rp), (35) 

Eqs. (j3"2")l and (|3"3")l can be written in the simpler form: 



2 Y^^(q 2j +i)v^ a2 P cos7rp 
j J p^O F J 



4 v-^ e|,(g 2 j+i)cos7rZ 



a 2l (E}{k 2l ) - hW) + —J2 

j 

y-\ a 2p f< ;t 7; 



7T- ^ 2/ - 27 - 1 

3 J 



p^O 



2p -2j-l 



After the substitution 

a 2p cos np = a 2p 

Eq. (|3"T)) becomes 

a 2l (E}{k 2l ) - H 2 uj 2 ) 



J_Y^ £ 2 /(92j+i) a 2p 



7T 2 21 - 2j - 1 2p - 2j - 1 



0. (37) 



(38) 



0, (39) 



where j and p are integers. We now take into account 
that, at any integer /, 



±y l — 

it 2 Z-j (21 -2j-l 



1. 



(40) 
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where j are integers (j = 0, ±1, ±2, . . .). At large N, 
the function el l (q 2 j+i) varies slightly, if I is changed by 
1. The main contribution to sum (|40j) is made by terms 
with 2j = 21 ± 1 and the nearest ones. Therefore, with a 
high accuracy, we have 



i 2 (01\ - AV £ 2i(ff2j+l) 

1 > ~ rr 2 4- (2/ - 2j - l) 2 



zfez)- 



(41) 



Let us separate the term with p — I in the sum in ([39} . 
Then, with regard for (l4"T} . we obtain the infinite system 
of linear homogeneous equations for the coefficients a 2p : 



4 \- 4l(l2j+l) 



IT 2 ^ 21 

j 



2.j 



P#0,i 



G2p 



2p - 2j - 1 



= 0. (42) 



The system has a solution, if its determinant is equal 
to zero. This condition yields the system of equations 
for ui 2 . Since ui 2 stands only at the diagonal elements 
of the matrix, every k 2 i corresponds to a single solution 
uj 2 . In other words, the number of frequencies lu 2 is equal 
to that of wave vectors k 2 i, and the dependence Lo 2 (k 2 i) 
sets a dispersion curve. If several collections of solutions 
{a 2p } exist, then we may have several dispersion curves. 

By changing the sign of / and those of (2j + 1) and p, 
we obtain from (|4"2"} : 



a- 2l (E 2 (k 2l ) - H 2 uj 2 ) 

IT 2 



21 - 2j 



i £ 

P#0,Z 



-2p 



2p -2j-l 



0. (43) 



Eqs. (j42|) . (|43[) imply that a 2p must satisfy the relation 

a_ 2p = za 2p , z = -1, (44) 



a-2p = za 2p , 



1 



(45) 



for all p. These arc two different classes of solutions. With 
regard for (44} and (|4"5j). Eq. ((42j) is reduced to 



«2i(^ 2 (fc 2 + eb(20 



p#Z,p>0 



1 



2p - 2 j - 1 2p + 2j + 1 



e 2i(g2j + l) 

2? - 2j - 1 
0, (46) 



;(20 



4z 



tt 2 ^ 4/ 2 - {2j + 1) 



(47) 



where l,p — 1, 2, 3, . . ., j — 0, ±1, ±2, . . .. By excluding 
the term with p = I from the sum in Q39p , the sum con- 
tains only the terms with alternating denominators. This 
favors the smallness of the sum, since 



E 



i 



(2Z-2j-l)(2p-2j-l) 



at p ^ I (48) 



for integer p and I. These properties allows us to find 
the following method of solution of Eq. (|46|) . We as- 
sume that the dispersion law h 2 uj 2 (k 2 i ) is known. Let 
the frequency correspond to the wave vector k 2 i : uj 2 = 
u 2 (k 2 i ) = uj 2 (21 ). From the system of equations (|46j) . 
we find a 2 i^ 2 i as functions of a 2 i and lu 2 and then sub- 
stitute these a 2 i in Eq. (|46|) with I = Iq; in it, a 2 i is 
canceled, and we obtain the equation for the frequency 
w 2 (2Z ). Eq. (46} yields 



a 2 i#2Z (2Z ) - „ 2 ,, -y— 2 



,4 2 ,(2Zo)^ 2 (2Z ) 



(49) 



E 2 (k 2l )+e 2 (2l)-h 2 uj 2 (2lo)' 
The quantity A 2 i(2l ) is sought by perturbation theory: 
A 2l {2l ) = A° 2l (2l ) + SA 2l (2l ), (50) 



i(?2i+l) 



a 2Z 



h 2 u; 2 (2l ) 2/-2j-l 



1 



2Z - 2j - 1 2Z + 2j + 1 



(51) 



P#±Z,±Z 



SA 2l (2l ) 



E 



e 2z(<?2j+l) 



^ 2 (fc2 P ) + e 2 (2p)-ft 2 W 2 (2^) 



21 - 2j - 1 



1 



2p - 2j - 1 2p + 2j + 1 



(52) 



The idea consists in the following. Since the frequency 
corresponds to fc 2 z , we consider the harmonic 21q to be 
principal and separate it in the sum in (|46l) , by assuming 
that the contribution of the other harmonics (with p =/= Iq) 
to the sum is small. In this case, A 2 i can be represented in 
the form ([50]) as the sum A^ ([51]) and a small correction 
SA 2 i (|52|) . This is true if, for all I and Iq, 



\SA 2l (2l )\ « |^(2Z )|- (53) 
For I — Iq, (|46|) yields the equation for the frequency: 
h 2 u 2 (2l ) = E 2 (k 2la ) + e 2 (2l ) + SE 2 (2l ) 1 (54) 



SE 2 (2l ) = - 

7T Z 



p>0,p^Z o 

E ■ 

1 



K 2I 



+ 1, 



2/n 



2j - 1 a 2to 
z 

_2p-2j - 1 2p + 2j + 1 
If ([53]) is valid, we have 

^ h 2 (2Zo) « -^ 2 w 2 (2z- )M 2;o (2/ )/a2Z . 



(55) 



(56) 



According to (52} , <5A 2;o (2z' ) - (fi 2 w 2 (2fo)) _1 . This re- 
sult and the structure SE 2 (2l ) (52} imply that <5 J B 2 (2/ ) 
should be normalized by e|, ((? 2 z ), rather than by 
h 2 uj 2 (2l ). Then 



<5^(2Z ) = g(2Z )el, (^ ), 



(57) 



5 



q(2l ) « - 



(k2i ) a 2 i 



(58) 



It is difficult to analyze analytically the equations ob- 
tained, but the numerical analysis is rather simple. For 
He A atoms, the numerical analysis indicates that con- 
dition ([53| is satisfied if Uoa/R < IK at z = — 1 and 
Uoa/R < 0.1 K at z = 1, which coincides with the Bo- 
golyubov criterion Uoa/R <C fi 2 /(2mo 2 ) (in the ID 
case), because a s» 2 A and h 2 /(2ma 2 ) fa 2 K for i?e 4 
atoms. Here, Uq is the height of the potential barrier, a 
is the width, and R is the mean interatomic distance. In 
the numerical analysis, we use a simple potential 



U(x) 



U o >0, r<a 
0, r > a. 



(59) 



Let us consider Eq. If we represent e 2 (|47|l as 

^(2M=4o(^ )C(Zo), (60) 
we find numerically that 

c(i )^±^(ty, a ~i, (6i) 



far from the zeros of the function e 2 io (fc2/ ). Near these 
zeros, |C(/o)| tends to 1; at a zero, |C(/o)| — > oo. But, 
the inequality |e 2 (2Z )| < E 2 (k 2to )/N holds at all points. 
Therefore, the quantity e 2 can be neglected in ((54)) (in 
other formulas, we retained it in numerical calculations). 
Then Eq. (|54"|) can be written in the form 



h 2 uj 2 (2l Q ) 



h 2 k 2 

2m 



+ (2 + q(2l ))n v(k) 



2 1.2 



2m 



(62) 



with q (|58|) . The numerical analysis indicates that 
\6A 2 i (2l )/a 2 i \ < 1 and \q(2l )\ < 1 for z = 1 and 
z = — 1. Finally, we obtain 



(63) 



It is the famous Bogolyubov formula (|2~4")l . We see that 
its derivation under the consideration of boundaries is 
much more complicated than that under periodic BCs. 

In the numerical analysis, we used h 2 ui 2 (|62[) with q 
as a free parameter. We varied it from -10 to 10 and 
compared with the theoretical q (fS"5|). The solution is 
such q, for which the theoretical q (|58")l is equal to a free 
one. We obtained q = ± 0.001 for N > 100 and for the 
maximum values of I, j,p > 100. If N, l,j,p increase, the 
results do not change noticeably. On the whole, q is very 
reliably determined. 



B. Second solution. 

While obtaining the first solution, we expanded the 
odd harmonics arising in Eq. (|29l) in the even ones in 



order to reduce all terms to the independent basis func- 
tions described by even harmonics. Such method leads to 
a definite combination of harmonics and to a certain dis- 
persion law. But we can obtain another combination of 
harmonics, which is also a solution. To make it, it is nec- 
essary that all harmonics arising in the equation be in the 
starting expansion of the functions nix) and s{x). Then 
the mutual compensation occurs for other collection of 
harmonics, and we may obtain a new dispersion law. The 
physical motivation is simple. The system admits stand- 
ing waves with k — (2j + l)ir/ L and k = 2jir/L. They are 
independent physically, and it is possible to attain their 
mathematical independence. For this purpose, we repre- 
sent h(x) and s(x) as the sum of two functions, where the 
first function is expanded in even harmonics e lk2lX and 
the second function is expanded in odd ones e lk2j+1 ' x : 



n{x) = 



a 2 ie 



t21x/L 



+ Z^ a 2j+ie 



iTr(2j+l)x/L 



3 



(64) 



s(x) = V b 2l e^' L + V b 2j+1 e i ^+ i y^ L , (65) 



3 



here I and j run all integers. The both collections of ex- 
ponents e lk2lX and e lk2j+lX are complete sets of orthogo- 
nal basis functions. We substitute these formulas in Eqs. 
(HU), PD|. ppj) . The integral in flU gives the sum of 
even and odd harmonics. We do not expand them in any 
other harmonics but gather them with the same harmon- 
ics from other terms. Gathering all coefficients of each 
of the functions e lk2lX and e lk 2j+i x ^ we equate the sums 
obtained to zero and arrive at the system of equations 
for the coefficients a n and b n . Relation (fj"9|) yields 



0, b 



no 



1,2 u ' n - 



(66) 



In view of (l66l). we obtain from (1201): 



1/(0)^ sin [n(2j + l)/2] 
ub = > a 2 j+i — — , (67) 

7]- £ 1 



2j + 1 



a 2 i(Ej(k 2l ) - h 2 co 2 ) 

2 2 fi. a 2 i+ i sin [tt(2j + l-2/)/2] 



(68) 



2j + l- 21 



0, 



a 2j+ i(E 2 (k 2j+1 ) - H 2 uj 2 ) 



(69) 



2 2 ^ ag; sin [tt(2Z - 2j - l)/2] 

where i, j = 0, ±1, ±2, . . ., but / ^ 0. Since 
sin [tt(2Z - 2j - l)/2] = - cos (ttZ) sin [7r(j + 1/2)], (70) 

it is convenient to make changes 

a 2i cos (irl) = a 2U a 2j+1 sin [ir(j + 1/2)] = a 2j+1 . (71) 
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Then Eqs. (|68l) and (|69j) are simplified: 



&2i(E 2 f (k 2l )-hW) + -e 2 2l (k 2l ) 



2j + l- 21 



= 0, (72) 



a 2 j+i{E 2 (k 2 j+i) - h 2 u> 2 ) 

The symmetry of the equations is such that 

a_22 = -za 2U a_2j_i = 2a 2i+1 , (74) 

where 2 = 1 (the solution A) or z = — 1 (the solution B). 
In view of this, we write (|72p and ([73]) as 



a 2l (E)(k 2l )-h 2 u?) 
2 

= 4(fe;)y]a2j+i 



(75) 



(f + z)2/+(l-z)(2j + l) 



(2j + l) 2 - 4Z 2 



a 2]+1 {E 2 f {k 2 . ]+1 ) - h 2 uj 2 ) 



(76) 



2 2 ^ (l + z)2/ + (l-z)(2j + l) 
- --%•+! (W) 2,021 (2j + l)*-4l* ' 

where l,j = 0,1,2,3,..., but I 7^ 0. Substituting a 2 j+i 
(|76|) in the equation for a 2 z, we obtain the equation only 
for the even coefficients: 

a 2l {E%k 2l ) + e)(2l)-hW) 

4 2 /1 \ - e lj+l( fc 2j + l) 
" ^ £ 2 ; (^)^ i? 2 (fc2j + i) _ ft 2 c ,2 

(l + z)2Z+(l-z)(2j + l) 

X (2j + l) 2 - 

v x (l + z)2/ 2 + (l-z)(2j + l) 
x 2, fl 2i 2 , ^„ 772 , (77) 



(2j + 1)2 - 4*1 



^ £ 2( 



(ME 



t 2j + l 



(<?2j+l 



E}(k 2j+1 ) - 

[(l + z)2l+(l-z)(2j + l)} 2 
[(2j + l)2-4F]2 



(78) 



where l,l 2 ,j = 0,1,2,3,... and l,l 2 ^ 0. In (77]), the 
term with l 2 = I is separated from the sum and is trans- 
ferred to the left-hand side. 

We solve Eq. ([77jl by perturbation theory described 
above. We separate the harmonic I — Iq as the basic one, 
so that uj 2 — u; 2 (2/ ). For I 7^ Iq, we have 



«2Z/2i (2?o) 



A 2/ .(2Z )Zl 2 u; 2 (2Zo) 



E 2 (k 2l )+e 2 (2l)-h 2 u 2 (2l )' 



A 2l (2l ) = A° 2l (2l Q ) + SA 2l (2l ) 



(79) 



(80) 



4 4(M 



E 



S2io e 2j+l(^2i+l) 



7T 2 h 2 oj 2 (2l ) ^ E 2 (k 2j+1 ) - h 2 uj 2 {2l ) 

3 •* 

[(l + z)2Z + (l-z)(2j + l)] 

[(2j + 1)2 - 4/2] 
[(l + *)2Z + (l-«)(2j + l)] 



[(2. ? + l)2-4/2] 



(81) 



5A 2l (2l ) » _4(fc 2i )E 



J 2j + l( fc 2j + l 



E 2 (k 2j+1 )~h 2 uj 2 (2l ) 



(1 + z)2l + (1 - z)(2j + 1) 



E 



(2j + l) 2 - 4P 

(l + z)2Z 2 + (l-z)(2j + l) 
(2j + l) 2 - 4/2 



^ 2 (2/o) 



E 2 (k 2l2 ) + e 2 (2l 2 )-h 2 Lu 2 (2l )-' 



(82) 



^ ^2, J = 0, 1, 2, 3, . . . and /, l 2 7^ 0. For the mode / = Z , 
we obtain the dispersion equation 

h 2 Lu 2 (2l ) = E 2 (k 2lQ ) + ej(2l ) + SE 2 {21 ), (83) 



-Ij+iihj+i) 



SE H 2l o) --A l0 (k2 la ) E E 2 {k 2 ^)-h 2 u> 2 (2l ) 
(l + z)2Z + (l-z)(2j + l) 



(2j + l) 2 - 4/2 
a 2h (l + z)2/ 2 + (l-z) (2.7 + I) 
a 2lo (2j + l) 2 - 4/ 2 2 



(84) 

L r - *'2 

Under condition (|53p. we have 

^(2/ ) w -^ 2 w 2 (2/o)M 2 / Q (2/o)/a 2;o . (85) 

Values of 8E 2 (21q) and e 2 (2Zo) are proportional to the 
quantity e 2 ; (& 2 ; ), rather than to H 2 lu 2 (21 ). For e 2 (2Z ) 
(f75|) . it is obvious. For SEj(2l ) ([55]). it is easy to verify. 
Therefore, we write 



4(2i ) = 4o(*ao)*i(2io), 



(86) 



^ 2 (2/ )-e 2 ;o (A: 2;o ) S2 (2/o). (87) 

Then the squared frequency ([83]) reads 

^ 2 (2/ ) = S 2 (fc 2/o )+e 2 ;o (fc 2/o )[ Sl (2/ ) + S2 (2/ )] 
/ \ 2 

/ H 2 k 2 , \ nv(k 2] ) h 2 k 2 , 
1 1 +g(2Z ) ^,(88) 



2m 



g(2/ ) = l + Sl (2/ ) + S2 (2/ ). 



(89) 



The numerical analysis of formulas (TT8"]) - (j5iJ[) indicates 
that, for He 4 atoms, the perturbation theory works well 
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for Uoa/R < 0.001 if. Moreover, the smallness of Uo 
is more significant than that of a/R (but this can be 
related to the fact that, as a/R decreases, it is necessary 
to sharply increase the maximum values of j, I2 in the 
sums, which causes a growth of the error of the method). 

We considered the parameter q in ([88]) to be free and 
varied it in wide limits. Then we calculated the theoret- 
ical q (|89| . It turns out that, for lo > N, the theoretical 
q is equal to 1 with high accuracy (|si|, | S2 ] < 0.01) for 
any bare q (as for the fictitious solutions, see below). As 
lo decreases to ~ 0.5 AT, the theoretical q starts to de- 
pend on the bare one. For l < 0.27V, this dependence 
becomes strong. The solution is such q, for which the 
theoretical q coincides with the bare one. Its dependence 
on l and N for a/R — 1 is given in Fig. 1 (for a/R <C 1, 
the solutions are similar). As is seen, q — 1 for large lo 
and various N. For small Iq, q depends strongly on N. 
For lo ~ 1, q decreases, as A~ increases, and attains a 
minimum at A^ ss 500. As A^ increases further, the value 
of q increases, by attaining the maximum for A~ 3000. 
With the further increase of N, q makes small oscillations 
near q = 1. For N > 10 4 , the equality q = 1 holds with 
high accuracy for all Iq, i.e., we have the solution 



h 2 u 2 {k) 



E 2 Jk) 



h 2 k 2 

2m 



nv{k) H 2 k 2 



(90) 



It differs from the Bogolyubov formula by the factor 1/2 
of the potential. Though we have found the solution 
only for k = ir2l/L, it is clear from the symmetry of 
Eqs. (|75|) and ([76")) that the solution is the same also for 
k = tt(2Z + 1)/Z. 




10 15 20 25 



Fig. 1: Values of q(l ) for the solution 2 (Eqs. (JSSj, Hi with 
k — ty21o/L, z — —1) for the system of N atoms He 4 with 
the interatomic potential (|59[l; a/R = 1 and Uo = 0.0033 if. 
Circles show the curve for N — 100, triangles — N = 500, 
crosses — N = 3000, rhombs — N = 8700, and stars - 
N = 16000. 

We note that the same dispersion law (l9Tfl) was ob- 
tained with regard for boundaries in more exact ap- 



proaches [3, [lj| , where the total AT-paxticle wave func- 
tions of a system were calculated. In |3 , the general 
formula for the Bose liquid in a rectangular box of any 
dimensionality was obtained: 



'h 2 k 2 
2m 



nv{k) h 2 k 2 
~2~f m~' 



(91) 



Here, / is the number of noncyclic coordinates of the 
system. In [l8|, [T9[, only one dispersion law (|9"Tj) was 
deduced. For / = 1, it passes into (|9"U1) . Our approach 
and that in pjl, HH are completely different. But both 
yield law (I9U1) . Hence, it is clear that this law is not an 
incidental error and must correspond to the physics of a 
system. 

In addition, there exist a huge number of fictitious so- 
lutions for every with different q. For them, the de- 
nominators in JEU) and (1521 are close to zero: 



£5(ty + l)-fiV(2lo)«0 



E 2 f (k 2l2 



h 2 u 2 (2l Q ) « 



(92) 



(93) 



They are distorted solutions of the 



with n 2 uj 2 {2lo) 
form (|50"|) , but with different k. Such q satisfy (1581 and 
(|89| . but condition (|53|) is not satisfied for them, i.e., 
we cannot use (J58J and . If the denominator in ([9^)1 
becomes zero, this means that we must choose some other 
mode with 2lo = 2l 2 as the basic one, and we will obtain 
law (|M|) for it. The equality of (|9"2")l to zero means that 
the odd mode 2j + 1 is the basic one, and we must solve 
the equation for odd harmonics ([To]) , by substituting ([75)1 
in it. We obtain again (f90|) . 

Though the deviation of q from 1 for small Iq and N 
is not usually related to fictitious solutions, the fictitious 
solutions affect noticeably the true solutions for some N 
(TV=3000 and 8700 in Fig. 1, JV = 3506, and many oth- 
ers), by shifting them. This is observed for crosses in Fig. 
1, for which the nearest fictitious root q s» —7.2 from 
class (1921 . and for rhombs (q ~ 3 from class (|9"2")l ) and for 
N = 3506 (the curve is similar to that from crosses, but 
with q(lo = 1) ~ —0.07; here, the fictitious q ~ 2 is from 
class (|93"l) h For A^ > 10 4 , there exist also many narrow 
strips of values of N, for which the fictitious root is lo- 
cated closely to 1 and deviates the solution, q, from 1 for 
lo ~ 1. In this case, the curve E(k) (|9"01 bends at k — > 0. 
As a result, the critical velocity v c = min(E(k)/Hk), at 
which the superfluidity is destroyed, is set by the domain 
of k — > and is much less than the sound velocity. For 
Iq 1, no deviations of q from 1 were found. 

We studied Eq. (JTTJ) and system ([75]). (175)) within an- 
other method: we determined numerically the modulus 
of the matrix determinant, by setting 



frV(fc) 



h 2 k 2 

2m 



ni>(ko 



h 2 k 2 



(94) 



where A: is the same for all q e '^ and is equal to ko- The 
values of q e ^ , for which the determinant tends to zero, 
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It is seen from Fig. 3 that, for small k, the results are 
almost independent of fco and are close to the above ones 
determined by perturbation theory. As k increases, the 
accuracy becomes worse. The method is exact for A: in a 
vicinity of ko. In order to obtain exact q for large k, it 
is necessary to select a large fco ~ k. The advantages of 
this method are the simplicity, the absence of fictitious 
solutions (because zero denominators are absent), and 
the possibility to calculate the frequencies both even and 
odd harmonics at once, see Fig. 3 (whereas, by pertur- 
bation theory, we must consider even and odd harmonics 
separately). As a drawback, we mention the absence of 
analytic formulas for the dispersion law. 

As the dependence q(k) (Figs. 1, 3) has been found 
within two methods, it is reliable. 

In the same way, by calculating a determinant, we have 
solved the system of equations (|46|) . The solution is iden- 
tical to that found by perturbation theory: it is given by 
formula {52} with q = ± 0.001 (for z = ±1). 



C. False solutions. 

At the substitution of ^ and (J53J) in Eq. with 
potential (I25[) , we obtain for example, terms of the form 
a 2 iv{k 2i )e lk2lX and a 2]+ iv(k 2 j + i)e tk2 ^+ lX . The former 
can be presented as 

a 2l v{k 2l )e lk2lX = r 2 (k 2l )a 2l v(k 2l )e lk2lX (96) 
+ (1 - r 2 {k 2l ))a 2l u{k 2l )e^ x = r 2 {k 2l )a 2l v{k 2l )e^ x 

+ (1 - r M) a 2 Mk.) E " i'^f ^ 



Fig. 2: Logarithm of the modulus of the determinant of 
the matrix of coefficients of Eqs. (|75[) and (|76[) versus q e ^ 
from (|94[) for z = — 1 and fco = n/L (circles), fco = 2n / L 
(stars). The system contains TV = 100 atoms He 4 with the 
interatomic potential (|59p with a/R — 1 and Uo = 0.0033 K. 



Fig. 3: Solutions q(lo) for the system of equations (|75|l . 
([76} with z = — 1 and fc = tvIq/L for atoms He* for the 
same TV, a/R, and Uo, as in Fig. 2. Triangles show the 
curve for fco = n/L in (|94|l . squares — for fco = 2tv/L, and 
the circles taken from Fig. 1 show the curve calculated by 
perturbation theory. 



are solutions. The characteristic dependence of the deter- 
minant on q e M is shown in Fig. 2. The sharp minima of 
the determinant indicate the roots of the equation. By 
them, it is possible to reconstruct the required depen- 
dence q(k) (Fig. 3) for various fc, by using the formula 



h 2 k 2 a 

2m , 



K 2 k 2 

2m 



q(k)- 



nv{k ) h 2 k^ 
2 m~ 

-v(fc) h 2 k 2 
2 m~' 



by expanding the even exponential function e lk2lX in 
odd ones e lfe 2j+ia; Similarly, we may do with the term 
a 2]+ iv{k 2]+l y k2 ^ x : 

a 2] +iv{k 2j+1 )e ik2 '+ lX = ri(fc 2j+ i)a 2j - + i^(fc 2j+1 )e lfe ^+ 1 
+ (1 - ri{k 2j+1 ))a 2 . ]+ ii>(k 2j+1 ) x 
^ sin 7r(Z-j - 1/2) ikmx 
4- -1/2) 



(97) 



(95) 



Then, instead of f|68[) and f|69[) . we will obtain another 
equations, which yield a lot of other dispersion laws. 
However, they are nonphysical. The presence of many 
solutions is related to the fact that the frequency is de- 
termined from the system of equations, rather than from 
a single equation. When we expand some harmonic in 
other harmonics in each equation, this transformation is 
equivalent for the equation, but not for a matrix. In the 
last case, the harmonics are redistributed in a row, and 
the determinant of the matrix is changed. Hence, the 
frequencies are changed as well. Thus, by transforming 
relations (l96l) and (|97|) , we replace the initial system of 
equations under solution by another system in an arbi- 
trary way. This leads to false solutions, and we will not 
study them. 
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III. GENERAL PROPERTIES OF SOLUTIONS. 



Solution 1 is independent of the initial collection of 
basis functions: the same equations are obtained, if we 
expand n(x) and s(x) in odd exponents e 2 *+ lX on the 
interval [— L/2, L/2] or in sines sin (Iwx/L) on [0, L]. For 
solution 2, we have not studied other expansions. 

Let us consider the properties of solutions for h(x) and 
s(x). For z = 1, solution 1 corresponding to the 2Zo-mode 
takes the form 



n 2i a {x) = 2} a 2 ;(2? ) cos (irl) cos (2-kIx/L), 



(98) 



S2i (x) = b (2l ) 



2muj 
n 



^ 2 ^ 2 °^ cos (irl) cos (2nlx/ L), 

l 2! 

(99) 

i = 1, 2, 3, . . ., k p = irp/ L. The coefficients a 2 /^2i are de- 
termined from formulas (H^l) - ([51|) . (jf)3"j) . where we assume 
5A2i(2Zo) = and a,2i = no- In Figs. 4 and 5, we present 
the values of 02; and the function n,2i (x) (|98p for two Zo- 
For z = — 1, the solutions are analogous, but cos (2nlx/L) 
in and should be replaced by i sin (2-kIx/L). It 
is also necessary to set a2i — —iriQ. Then bo = 0, as 
follows from (|36|) . For z = — 1, the function n2i (ir) is 
equal to zero on the boundaries (n2i (x = ±L/2) = 0). 
For z — 1, it is different from zero. 

Solution 2 corresponding to the Zo-mode for z = — 1 
takes the form 

ftj (a:) = 2^ a 2i (Z ) cos (ttZ) cos (2ttZ.t/L) (100) 
l 

+ 2^2 SiW+iCo) sin [tt(j + 1/2)] cos [tt(2j + l)x/L], 



S; (x) = 6 (Z ) 



2mu; x - a 2 i(Zo) 



1,2 
2/ 



cos (ttZ) cos (2-kIx/L) 



(101) 



2mu> 
n 



E 



Q2j+l(^o) 
K 2j + 1 



sin [7r(j + 1/2)] cos [n(2j + l)x/L], 



where ^ = 1,2,3,... and j = 0, 1, 2,3,.. ., and the coeffi- 
cients ai are given by the formulas ([75]). ([75 ]) -([8T| ) . (1581 . 
and ([59^1 . There, we set 5 A = and &2; = no. Relations 
([g?]). CD and GH yield the real 6 . In Figs. 6 and 7, we 
show the values of a; and the function hi (x) (|100p . 

For z = 1, the solutions are the same, but cosines of 
the form cos (ai) in (|100|) and (|10ip should be replaced 
by i sin (ax) with the imaginary unity i. To get rid of 
the imaginary values, we set fi; = —mo. Since all terms 
in (llOOp and (|101|) are proportional to a; , the imaginary 
unity disappears everywhere. According to (j57J, (J7TJ) and 
(1741 . we have b = 0. 

The figures show that, for solutions 1 and 2, the wave 
packets are localized in the fc-space and are delocalized 




15 20 25 30 



I 




Fig. 4: Dependence a,2i(2lo)/no on I for solution 1 for z = 1 
and Zo = 3 (rhombs), Zo = 10 (circles). The calculations 
were performed for the system of N = 100 atoms He 4 with 
the interatomic potential (|59[) with a/_R = 0.1 and Uo = 
0.001 K. The dashed line corresponds to &21 = 0. For main 
harmonics (two points with lOOO^/no = 1: 1 — 3 and 
Z = 10), the value of &21 is larger by factor of 1000, then it 
is shown in the figure. 



Fig. 5: Dependence of the density fi2i (|98|) on the coordi- 
nate x for solution 1 for z — 1 and Zo = 3 (circles), Zo = 10 
(rhombs). The parameters are the same as in Fig. 4. 



in the cc-space. Due to the strong localization in the 
fc-space, the function h(x) is close to a monochromatic 
wave ± cos (2irlox/L) for solution 1 and to ± cos (irlox/L) 
for solution 2. It is of interest that solution 1 corre- 
sponds actually to only single harmonic with I = Iq (as 
for the cyclic BCs), the rest of the harmonics are strongly 
damped. Though the system is placed on the interval 
[—L/2, L/2], we show only the interval [— L/2, 0] in the 
figures, since the function on the interval [0, L/2] are the 
same: h(—x) = n(x). 
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Fig. 6: Dependence of ai(lo)/no on Z for solution 2 with 
z = -1 for Z = 6, q = 0.853 (rhombs) and Z = 20, g = 0.987 
(circles). Here, g are solutions of Eqs. (|88[) and ((89}, and 
the rest parameters are the same as in Fig. 4. 

Fig. 7: Dependence of the deviations of the density hi (| 100p 
on the coordinate x for solution 2 with z = — 1 for Zo = 6, 
g = 0.853 (circles) and Z = 20, g = 0.987 (rhombs). The 
other parameters are the same as in Fig. 4. 



The general form of solutions 1 and 2 is as follows: 

no(x,t) = ^2cihi{x)T l n {t), s (x,t) = y]disi(x)T l 8 (t), 
i i 

(102) 

where c/ and di are any numbers, I = 2, 4, 6, . . . for solu- 
tion 1 and I = 1, 2, 3, ... for solution 2. 



the interval a; £ [— L/2, L/2], which are obtained without 
regard for the sewing of ^(x,t) on the boundaries. As 
usual, the discreteness of k arise due to the sewing. Be- 
cause our k wp are already discrete, solutions 1 and 2 are 
valid not for all BCs, i.e., these solutions are not general. 
For the general solution, k wp must be continuous. For 
zero BCs, the general solution will be studied in the sub- 
sequent work. This solution satisfies more complicated 
equations. In the present work, we restrict ourselves by 
the simpler solutions 1 and 2, since they are easy to be 
studied. We can make it comprehensively within two dif- 
ferent methods, which allows us to test the validity of the 
main result: the existence of two dispersion branches. 

Solutions 1 and 2 hold in the case where the barrier 
height U s is finite, and the number TV of atoms is suffi- 
ciently large. Indeed, the value of ^(x, t) on the wall is 
nonzero at a finite barrier and is determined by the bar- 
rier height and parameters of an atom (the interatomic 
interaction gives only a correction). The quantity h(x) 
can be estimated, if we associate a variation in the density 
to a phonon, by equating the energy N(3x 2 of an elastic 
wave to the phonon energy fku. We obtain that the am- 
plitude of oscillations of the density n(x) is proportional 
to R/ V-N, if the wave corresponds to a phonon with high 
k ~ tv/R, and to n(x) ~ R/N, if the wave corresponds 
to a phonon with low k ~ tt/L. At sufficiently large 
N, the amplitude of oscillations is small, and we have 
\SR(x)\ < |*o (x,t) | on the wall (since h{x) w 2R 6R). 
In this case, the boundary conditions are satisfied due to 
^o(x, t), and values of SR(x), h{x) can be arbitrary. The 
analogous situation arises in the approach [l8| with TV- 
particle wave functions, where the zero BCs are satisfied 
due to 4'o, so that BCs are free for phonons. 



V. POINT INTERACTION. 

We note that, for a system with point interaction, the 
exact solutions under periodic and zero BCs were deter- 
mined by Lieb, Liniger (20| and Gaudin (2lj . and no in- 
fluence of boundaries on bulk properties was found. Let 
us try to understand it, using the harmonic analysis. The 
effect appears due to the different ways to combine the 
harmonics and their different "weight" . The last circum- 
stance is related to the difference of v(k) for different k. 
For point interaction, an interatomic potential has the 
form of the (^-function. The (5-potential is particular. It 
is "bad" from the viewpoint of smoothness, but its formal 
expansion in a Fourier series, 



U5(xi - x 2 ) 



U \ " r i7Tl(x 1 -X 2 )/L 

i 



(103) 



IV. BOUNDARY CONDITIONS. 

The above-obtained solutions 1 and 2 describe the 
wave packets centered at discrete k of the form k wp = 
Tvj/L or 2irj/L. They are the solutions for &(x, t) on 



(I are integers), corresponds [22j to the 6- function by its 
properties. In (|103[) . all v(k) are identical and equal U. 
Such a Fourier-expansion has no dispersion. Therefore, 
other combination of harmonics does not lead to a new 
laww(fe). Indeed, let us return to the Eq. (|2^|) . Relations 
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p6t. d3l[). and (|35[) yield 



2a 2 



dx'h(x')U(\x - x'\) = ^M^H/l 



i2irpx j L 



]2,P 



cos np cos nl 



(2l-2j 2 -l)(2p-2j 2 -l) 



(104) 



We take the absence of a dispersion into account, 
v{k2j 2 +i) = U = ^(0), and separate the term p = j 2 
in the sum. Then we have 

dx'h(x')U(\x - x'\) = ^e ,2l,1 / L 



2a 2 ;^(0) i27r i x /L 



E 



(21 - 2j 2 - l) 2 



(105) 



2a 2 iv(Q) p ^ ei27r „ , , 



E 



cos np cos ttI 



(2l-2j 2 -l)(2p-2j 2 -iy 



In view of (gOJ) and (05]), we obtain 

dx'n{x')U(\x - ar'|) = a 2l v{0)e l27rlx/L = v{Q)h{x). 

(106) 

The same result is obtained for periodic BCs for any po- 
tential. In other words, for a point interaction, the ape- 
riodic and periodic boundaries give the Bogolyubov dis- 
persion law. This corresponds to the results [2(| HH by 
which the energies of a cyclic system and a system with 
boundaries coincide (with a negligibly small difference). 
However, the interactions are nonpoint in the Nature, 
and our solutions hold for a nonpoint interaction. 



VI. CONCLUSION. 

We have shown that, under the consideration of bound- 
aries, the Gross-Pitaevskii equation has two solutions for 
the dispersion law of phonons: the Bogolyubov solution 
and a new one with a factor of 1/2 of the potential. Both 
dispersion laws correspond to the same ground state. Be- 
cause only a single solution is realized in the Nature, one 
of the dispersion branches must be unstable. In work 
[23j . the Gross-Pitaevskii equation was solved for a vor- 
tex ring, and two dispersion laws were also obtained. The 
upper branch turns out unstable. We have not studied 
the stability of the curves. Since the new curve corre- 
sponds to a lower energy, this curve is probably stable. 
It is not easy to clearly explain why the boundaries cause 



the appearance of a new mode. The reason is related to 
the topology, namely to the difference in the collections of 
eigenharmonics for closed and open systems (see [l8[ for 
details). Moreover, two curves should correspond to two 
different ways to the diagonalization of a Hamiltonian. 

The author is grateful to Yu. V. Shtanov for the discus- 
sion of the results. This work is partially supported by 
the Special Program of Fundamental Studies of the Divi- 
sion of Physics and Astronomy of the National Academy 
of Sciences of Ukraine. 



VII. APPENDIX. THE CONDENSATE IN A ID 
SYSTEM 



The following formula was obtained [15[ for Green's 
one-particle function in ID: 



G[x — x' '; i)|| x _ x /|_ 



-inoe 



-ykT\x — x' I / p s 



(107) 



(7 is the thermodynamic constant). For an infinite sys- 
tem, the largest \x — x'\ is equal to infinity. Therefore, at 
any finite temperature, the right-hand side of (|107[) tends 
to zero, which means the absence of both a long-range or- 
der and the condensate in ID. But, for a finite system, 
the maximum value of \x — x'\ equals L (system's size). 
Therefore, as T — > 0, the right-hand side of (|107p tends 
to —mo, i.e., the condensate is present. For a ID sys- 
tem, the number of phonons N p h/L ~ T. For the least T 
(at which the system contains a single phonon), we ob- 
tain T m i n ~ 1/N. In this case, \x — x'\ « L ~ N. Hence, 
jkT\x-x'\/p s ~ 1, saidG(x-x'-t) -in C,C - 1. Of 
course, the temperature is defined only for a large num- 
ber of quasiparticles. But it is clear that, for low-excited 
states, the condensate is not forbidden (analogously, we 
can verify that, in a 2D system, the condensate is per- 
mitted up to temperatures comparable with T ~ IK, 
when the superfluidity is broken). For the other proofs 
0, El, Q3 > the same conclusion and an analogous esti- 
mate for T are valid, since the integral J dkT/k 2 does not 
diverge for a finite system (it can be replaced by the sum, 
where the zero level k — is transferred to the conden- 
sate, and the next levels k = 2tt/L, 47r/L, . . . are discrete 
and give a finite number for a finite L). In this paper, we 
consider that T is sufficiently low, and the condensate is 
present. 

Bogolyubov formulas 1] imply that, in a finite ID sys- 
tem, J2 N k/N < (a/RjlnN < 1 at T = (a is the size 
fe>o 

of an atom), i.e., the condensate is permitted at T = 
too. 
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